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ABSTRACT 

We quantify the morphological evolution of z ~ 0 massive galaxies (M*/M 0 ~ 10 1L2±0 ' 3 ) from z ~ 3 in 
the 5 CANDELS fields. The progenitors are selected using abundance matching techniques to account for the 
mass growth. The morphologies of massive galaxies strongly evolve from z ~ 3. At z < 1, the population well 
matches the massive end of the Hubble sequence, with 30% of pure spheroids, 50% of galaxies with equally 
dominant disk and bulge components and 20% of disks. At z ~ 2 — 3 however, there is a majority of irregular 
systems (~ 60 — 70%) with still 30% of pure spheroids. 

We then analyze the stellar populations, SFRs, gas fractions and structural properties for the different mor¬ 
phologies independently. Our results suggest two distinct channels for the growth of bulges in massive galaxies. 

Around ~ 30 — 40% were already bulges at z ~ 2.5, with low average SFRs and gas-fractions (10 — 15%), 
high Sersic indices (n> 3 — 4) and small effective radii (R e ~ 1 kpc) pointing towards an even earlier formation 
through gas-rich mergers or violent disk instabilities. Between z ~ 2.5 and z ~ 0, they rapidly increase their 
size by a factor of ~ 4 — 5, become all passive and slightly increase their Sersic indices (n ~ 5) but their global 
morphology remains unaltered. The structural evolution is independent of the gas fractions, suggesting that it 
is driven by ex-situ events. 

The remaining 60% experience a gradual morphological transformation, from clumpy disks to more regular 
bulge+disks systems, essentially happening at z > 1. It results in the growth of a significant bulge component 
(n ~ 3) for 2/3 of the systems possibly through the migration of clumps while the remaining 1 /3 keeps a rather 
small bulge (n ~ 1.5 —2). The transition phase between disturbed and relaxed systems and the emergence 
of the bulge is correlated with a decrease of the star formation activity and the gas fractions suggesting a 
morphological quenching process as a plausible mechanism for the formation of these bulges (although the 
eventual impact of major mergers and a growing black hole in the bulge should also be considered). The 
growth of the effective radii scales roughly with H(z)~ i and it is therefore consistent with the expected growth 
of disks in galaxy haloes. 

Subject headings: galaxies:evolution, galaxies:high-redshift, galaxies:structure 


1. INTRODUCTION 

In the local Universe massive galaxies are characterized by 
having a dominant early-type, bulge-dominated morphology 
as well as old stellar populations. They are also confined 
to tight scaling r e lations, such as the m ass-size relation (e.g. 
IShen et akl 12003b iBernardi et al.1 1201 4b and the fundamental 
plane. Understanding the formation and subsequent mass as¬ 
sembly of such systems is however still debated in present- 
day cosmology and it is a key milestone towards reaching a 
complete view of structure formation and the interplay be¬ 
tween baryons and their dark-matter hosts. In particular, the 
actual role played by mergers as compared to in-situ processes 
in shaping spheroids is still unclear, and state-of-the-art semi- 
analytic models of galaxy formation offer sometimes qui te 
different views (e.g jGonzalez et aLll2009HLapi et al1l2011l) 

Following the scaling relations of these massive objects 
across cosmic time is a natural way to better understand how 
the relations actually emerged and the role played at different 
cosmic epochs by the different physical mechanisms. As a 
matter of fact, many works in the last ten years have focused 


their attention on the evolution of the mass-size relation for a 
selection of massive galaxies (Zog(M*/M 0 ) > 10.5) finding 
an apparent increase of the z ero point of the relation by a fac- 
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changes in the slope (e.g . iNewman et al.ll20l2h or the scatter 
(Ivan der Wei et alJl20141) . 

Properly interpreting these redshift-dependent evolutionary 
trends as a progenitor-descendant relation remains still elu¬ 
sive given the continuo us mass build-up (e.g. iMuzzin et al.1 
12013: llibert et all 120131). the morphological transformations 


(e.g. Buitrago et al. 2011) and the evolution of the stellar pop¬ 
ulations (e.g IPeng et al.ll2010t) which make galaxies coming 
in and out of any given selection (progenitor bias effect, e.g., 
INewman et alJl2012l ICarollo et alJl2013l IShankar et al.ll2015l 

and references therein). As a matter of fact, a selection done 
at fixed stellar mass as usually done in several works, will nec¬ 
essarily be contaminated by galaxies which grow in mass that 
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will enter any given stellar mass bin. The level of contamina¬ 
tion depends on the stellar mass selection. For massive galax¬ 
ies (~ 10 n ), the fraction of galaxies in the lowest redshift bin 
which are actually descendants of the galaxies at higher red- 
shift (z ~ 2) is less than 20% (iMundv et al.ll2015l) . Therefore 
establishing evolutionary links is not straightforward at all. 
The situation is even worse when passive galaxies are to be 
considered since quenching and mass growth both contribute 
to this progenitor bias effect. 

One popular solution has been to study the evolution 
of the number density of these compact objects (e.g . 


Truiillo et al.ll201 ll: ICassata et al.ll20131: IPoggianti et al.ll2013b 


der Wei et al.l2014tlDamianov et alj2015l). but the results 

are not always in agreement, specially at low redshifts where 
HST surveys probe a small area and also because of the mul¬ 
tiple available definitions of compact galaxies. As a matter 
of fact, some works do select only the most massive galaxies 
(> 10e.g. iTruiillo et alJl201 fl) while others select 
a wid er stellar mass bin (> 10.5 M*/Mq, e.g IPoggianti et al] 
12013 0, On the other hand, there are papers using a fixed 
size threshold (a_ straight line in the mass-size plane, e.g. 
ICarollo et al.ll2013h while others prefer a s election according 
to the slope of the mass-size relation (e.g Ivan der Wei etalJ 
120141: ICassata et al. 11201 iBarro et ahll2oi3h . On top of this, 
other parameters that could bias the results are environment 
and a l so the passive selection criteria (e.g iPoggianti et alJ 
1.2013; IValentinuzzi et alJ l2010h . As a result, several au- 
thors do find a steep decrease of their abundances (e.g. 
ICassata et all2013Uvan der Wei et al .1120141) while other s tend 

to find a rather constant number (e.g. IPoggianti et al.1120131 : 

ICarollo et al. 11201 3l iDamianov et al.ll2015h . 

Another option has been propose d based on selecting galax¬ 
ies at fixed number density (e.g. Ivan Dokkum et al] 120101 : 
iPatel et all 1201 3l and references therein), i.e. assuming the 
rank order is preserved at all epochs. This approach also im¬ 
plies some important assumptions such as neglecting the role 
of mergers or the scatter in the mass accretion histories and 
it is faced to the known uncertainties in the evolution of the 
massive end of the mass function (e.g. Shan kar et alJ 120141 : 
iBernardi et allI 20 T 3 I 1 and the continuous quen ching happen¬ 
ing at all cosmic epochs (e.g. tPeng et al. 2010). Nevertheless, 
the latter approach can still provide some broad insights into 
the expected, average mass-growth of galaxies, thus allow¬ 
ing for a basic technique to observationally relate progenitors 
and descendants. Moreover, the methodology has now been 
improved by including corrections to the redshift-dependent 
number densities to account for mergers dPanovich et alJ 
l2014t iMarchesini et alJl20T4l : iBehroozi et al JI2013l) based on 
abundance matching. It was also recently stressed that dif¬ 
ferences in the stellar mass function lead to consistent results 
for the mass growth within ~ 0.25 dex, at least for galaxies 
with log(M*/M&) <11 (e.g.. iPapovich et afll2014l) . Glob¬ 
ally, these empirical studies based on number conservation 
procedures tend to agree on a significant structural evolu¬ 
tion, and confirm an important size growth of the average 
population. The growth seems to be driven by the addition 
of material in the outskirts of the galaxies (e.g. iPatel et aT] 
12013 0 in what has been called an inside-out growth and in¬ 
terpreted as a minor merger driven growth through the tidal 
disrup tion of small companions falling into the central galaxy 
(e.g. iNaab et al.l 120091: lOser et all 120101: iNinoti et alJl2012l : 
Shan kar et alJl2013l and references therein). 

Despite the outstanding efforts made so far, it is still chal¬ 
lenging to properly follow the evolutionary tracks of es¬ 


pecially the most massive galaxies. Number conservation- 
based approaches map today’s red and dead early-type sys¬ 
tems to progenitors presenting a variety of morphologies 
and star-formation activ it ies (e.g. ly tmdgrWdLgtid]|2008|; 
van Dokkum et al.l 1201 ll : iBruce et alJ 1201 2t IPapovich et al] 
2014t iMei et al.l l2014l> . It is thus very difficult to interpret 
the evolution of the overall population as a unique physical 
mechanism since multiple processes, such as morphological 
transformations or quenching can clearly contribute to move 
galaxies in the mass-size plane from one redshift bin to an¬ 
other. 

By simultaneously following the evolution of the star for¬ 
mation activity (quenching), morphological transformations, 
and mass build-up along the progenitor tracks identified 
through number-conserving techniques, one should be able to 
ideally separate the different contributors to the average struc¬ 
tural evolution. 

All previous works lack however of a precise quantification 
of how the morphologies change and evolve among the pro¬ 
genitors of massive galaxies, mainly. The most significant ef¬ 
fort has been probably done bv lBruce et al.l ( 2012 1 who made 
bulge-to-disk decompositions but just on one CANDELS 
field (UDS) and without quantification of irregularities^ An¬ 
other noticeable effort has been carried out bv lMortlock et alJ 
(120131) . who in the same UDS CANDELS field visually classi¬ 
fied galaxies into discs, ellipticals and peculiar systems. They 
found significant evolution in the fractions of galaxies at a 
given visual classification as a function of redshift, though 
they did not attempt to trace evolutionary tracks among galax¬ 
ies at different epochs. 

This is therefore the main new ingredient which moti¬ 
vates the present paper, in which we bring into the puzzle 
of massive galaxy formation detailed visual like morpholo¬ 
gies for a large sample of galaxies from all the five CAN- 
DELS fields. Combined with accurate stellar-masses and 
rest-frame colors as well as optical rest-frame structural pa¬ 
rameters from the 3D-HS T (iBrammer et al .112012| ) and CAN- 
DELS (iKoekemoer et al.l 1201 It Grogin et al.l 1201 ll) surveys, 
we revisit the evolutionary tracks of massive galaxies from 
z ~ 3. 

The paper proceeds as follows. In section [2] we describe 
the dataset used as well as the main physical parameters de¬ 
rived (morphologies, structural parameters, SFRs etc..). In §[3| 
we describe the procedure to select the main progenitors and 
from §0to §[6]we describe the main results, namely the evolu¬ 
tion of the morphologies, structures and star-formation prop¬ 
erties. These results are discussed in § [7] and we provide a 
summary in § [8] 

Throughout the paper, we adopt a flat cosmology with 
VLm = 0.3, S2a = 0.7 and Hq = 70 km.s~ l .Mpc~ l and we use 
magnitudes in the AB system. 

2. DATASET 
2.1. Parent sample 

We select all galaxies in the F160W filter with 
F160W<24.5 mag (AB) in the 5 CANDELS fields (UDS, 
COSMOS, EGS, GOODS-S, GOODS-N). Our starting-point 
catalo gs are the CAN DELS public photome tric catalogs for 
UDS (iGalametz et alJlMl and GOODS-S dGuo et alJl2013l) 
and preliminary CANDELS catalogs were used for COS¬ 
MOS, EGS and GOODS-N (private communication). The 
magnitude cut is required to ensure reliable visual mor¬ 
phologies (iKartaltepe et alJ l2014lb and structural parame- 
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ters ( Ivan der Wei et al. 112012l) which are two key ingredients 
for the analysis present ed in this work. As discussed in 
Ivan der Wei et all (i2014b . the magnitude cut results in a rea¬ 
sonable mass completeness of /og(M*/Mg) ~ 10 at z ~ 3 
which is well beyond the mass limit required to follow the 
progenitors of massive galaxies as discussed in the following. 
Our results should thus not be significantly affected by incom¬ 
pleteness. The median redshift of the sample is z ~ 1.25. 

2.2. Morphologies 

Visual-like morphologies are taken form the deep-learning 
morphology catalog described in Huertas-Company et al. 
(2015a) on the 5 CANDELS fields UDS, COSMOS, GOODS- 
N, GOODS-S and EGS. The classification m imics the CAN- 
DELS visual classification scheme from Kartalten e et alJ 
(120141) which is currently available in only one field. Mor¬ 
phologies are estimated using ConvNets, a specific artificial 
neural network t opology that is inspired by the biological vi¬ 
sual cortex (e.g. iFukushimall 19801) which is by far the most 
powerful image classifier up to date. When used for image 
recognition, convolutional neural networks consist of multi¬ 
ple layers of small neuron collections which look at small 
portions of the input image, also called receptive fields. The 
results of these collections are then tiled so that they overlap 
to obtain a better representation of the original image; this is 
repeated for every such layer. More details can be found in 
Huertas-Company et al. (2015a). 

The algorithm is trained on GOODS-S for which visual 
classifications are publicly available and then applied to 
the other 4 fields. Following the CANDELS classification 
scheme, we associate to each galaxy 5 numbers - f sp h , fdisk , 
fiff , fps, func - measuring the frequency at which hypothet¬ 
ical classifiers would have flagged the galaxy as having a 
spheroid, having a disk, presenting an irregularity, being 
compact or point source and being unclassifiable/unclear. As 
shown in Huertas-Company et al. (2015a), ConvNets are 
able to predict the fractions of votes given a galaxy image 
with a bias close to zero and ~ 10% — 15% scatter. The 
fraction of miss-classifications is less than 1%. We refer the 
reader to the aforementioned work for more details on how 
the morphologies are estimated. The important information 
to keep in mind for this work is that the classification is 
very close to a purely visual classification. We use only a 
classification in the H band (F160W) since the differences in 
the derived (broad) morphologies when using other filters are 
very small as shown in lKartaltepe et al.l 120141) . 

We are interested in distinguishing bulge and disk growth 
so we use the 5 morphology estimators to define 5 main mor¬ 
phological classes as follows: 

• pure bulges [SPH]: f sp h > 2/3 AND fdisk < 2/3 AND 
firr< 1/10 

• pure disks [DISK]: f sp h < 2/3 AND fdisk > 2/3 AND 
firr< 1/10 

• disk+sph [DISKSPH]: f sp h >2/3 AND fdisk >2/3 
AND f irr < 1/10 

• irregular disks [DISKIRR]: fdi S k >2/3 AND f sp h < 
2/3 AND f irr > 1/10 

• irregulars/mergers[IRR]: fdisk <2/3 AND f sp h <2/3 
AND f irr > 1/10 


The classification accounts for the presence or not of a 
disk/bulge component as well as asymmetries in the light 
profile. The thresholds used are somehow arbitrary but 
have been calibrated through visual inspection to make 
sure that they result in different morphological classes (see 
also lKartaltepe et al.l2014l) . We emphasize that slight changes 
on the thresholds used to define the classes do not affect the 
main results of the paper. Figure Q] shows some examples of 
the morphological classes defined that way. The SPH class 
contains galaxies fully dominated by the bulge component 
with little or no disk at all. The DISK class is made of galax¬ 
ies in which the disk component dominates over the bulge. 
Between both classes, lies the DISKSPH class in which we 
put galaxies with no clear dominant component. Then we dis¬ 
tinguish 2 types of irregulars: DISKIRR, i.e. disk dominated 
galaxies with some asymmetric features and IRR, which are 
irregular galaxies with no clear dominant disk component (in¬ 
cluding mergers). These two last categories do contain all 
the variety of irregular systems usually observed in the high 
redshift universe (e.g. clumpy, chain, taphole etc..). The sepa¬ 
ration between the last 2 classes is however challenging (even 
for the human eye), since a diffuse light component can be 
easily interpreted as a disk. Therefore, even though we will 
consider the 2 classes separately in most of the plots, the 
reader should keep in mind that there can be significant over¬ 
lap. For the galaxies selected in this work (see section [3]), 
> 95% of the population fits in one of the 5 defined classes. 
The remaining ~ 5% contains either galaxies with rather high 
irregular, spheroid and disk morphologies simultaneously or 
unclassifiable objects. 

2.3. Stellar masses and star formation rates 

Photometric redshifts, stellar masses and star formation 
rates (SFRs) are estimated from SEP model i ng as de¬ 
scribed in previous work s by IWuvts et~aTI (1201 lL 1201 2l) and 
iBarro et all (1201 3l 120141) . We describe here the basic pro¬ 
cedure and refer the reader to these works for more details. 
Photometric redshifts are estimated from a variety of dif¬ 
ferent codes available in the literature which are then com¬ 
bined to improve the individual performance. The tech¬ 
nique is fully described in iDahlen et al.l (120131) . Based on 
the best available redshifts (spectroscopic or photometric) we 
then estimate stell ar masses and UV-based SFRs using FAST 
(iKriek et al.ll20 09j) assuming iBruzual & Charlotl (120031) mod¬ 
els, a[Chabrier (2003) IMF, solar met allicity, exponentially 
declining star formation histories, and a lCalzetti et alJ (l2000h 
extinction law. Rest-frame magnitudes (U,V,J) based on the 
best-fit redshifts and stellar templates were computed using 
EAZY (iBrammer et al j|2008i) . 

The final SFR used in this work combines IR-based and 
UV-based (from SED fitting) SFRs as described in Bar ro et all 
(1201 lblal 120141) . The method essentially relies on IR-based 
SFR estimates for galaxies detected at mid- to far-IR wave¬ 
lengths, and SED-modeled SFRs for the rest. For Un¬ 
detected galaxies the total SFRs, SFRIR+UV, are computed 
from a combination of IR and rest-frame UV luminosities 
(uncorrected for extinction) following iKennicutl (119981) and 
iBell et alJ (120051) : 

SFRuv+ir = 1-09 x 10 w (L IR + 3.3L2&oo)[MQ.yr~ 1 ] (1) 

2.4. Structural properties 

Structural properties (effective radii, Sersic indices and 
axis ratios) are taken from the public catalog released 
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Fig. 1.— Example stamps of the different morphological types defined in this work from the COSMOS field. From top to bottom, spheroids, disks, 
disk+spheroids, asymmetric disks and irregulars. Coordinates are indicated in each postage stamp. 


in Ivan der Wei et afl d2012l) . Single Sersic 2D fits were per¬ 
formed to galaxies in CANDELS in the three infrared fil¬ 
ters (fl05,fl25,fl60) using galfit (iPeng et al.l 12002 ). The 
typical uncertainty on the parameters is less than 20% for 
th e magnitude cut applied i n this work as clearl y shown 
in Ivan der Wei et al.l (l2012l) . iBernardi et al.l (1201 4l) showed 
however that the total light profiles and sizes of massive galax¬ 
ies at z ~ 0 can significantly be affected by the background 
estimates. We do not expect a major impact of this effect in 
our sample at higher redshift though, where the contribution 
of the diffuse ligh t aroun d massive galaxies is less important. 
Ivan der Wei et al.l (120141) applied some corrections to the ef¬ 
fective radii of passive and star-forming galaxies to measure 
them in a unique rest-frame band of 5000. Given that the 
corrections are very small and have little effect on the final 
measured structural evolution as discussed in the aforemen¬ 
tioned work, we use here for simplicity the closest filter to the 
optical rest-frame band as done bv lNewman et ali (120121) . 

3. SELECTING THE PROGENITORS OF MASSIVE 
GALAXIES 

One key issue when one tries to infer the evolution of in¬ 
dividual galaxies is to actually link progenitors and descen¬ 
dants without being s trongly affected by prog e nitor bias (e.g. 
Carollo et al.l 1201ISonnenfeld et al.l l2014t IShankar et al.l 
20151) . The stellar mass function (SMF) is known to s ignifi- 
cantly evolve from z ~ 3 —4 (e. g. iPerez-Gonzalez et all 2008 ; 
iMuzzin et afl 120131: lllbert et al. 1201 3l) so a selection at fixed 
stellar mass will clearly be affected by new galaxies kicking 
in at lower redshifts as widely discussed in the recent litera¬ 


ture. Also, a selection of only passive galaxies will be affected 
by the continuous quenching at all cosmic epochs. An alterna¬ 
tive which is rapidly becoming very popular in the community 
is a se lection at fixed number density fe.g. Iva n Dokkum et al 


201 ftfBezanson~et al.ll201 ll: IConselice et al. 2013: iPatel et al 


2013). At first level, this selection assumes that the rank¬ 
ing of galaxies is preserved at all redshifts and therefore 
deliberately ignores the impact of mergers and the scat- 
ter in the mass accretion histories dBehroozi et al.l 12013 !) 
which can lead to errors in the stellar mass grow th of 
d(logM*)/dz ~ 0.16 dex (see also iLeia et aD 1201 3l for a 
comparison with SAM predictions leading to similar con¬ 
clusions). To overcome this issue, iBehroozi et al] d2013l) 
used abundance matching techniques to track the evolution 
of galaxies within their dark-matter haloes and apply a cor¬ 
rection to this simple assumption. The model therefore ac¬ 
counts for number density evolution and is the one adopted 
in this work. Figure [2] shows the stellar mass growth track 
for the progenitors of ~ 10 11 2 galaxies from j ~ 4 from 
the Behroo zi et alJ (1201 3l) model, assumin g the stellar mass 
functio ns (SMFs) of Baldrv et alJ (1200 81): iMoustakas et afl 
(1201 3l): fP erez-G onzalez ei al.l ( 20081) : iMortlock et al.l (1201 ll) : 
Marchesini et afl (120091 l20Toi) ~ As recently shown by 
Panov ich et afl ( 20141), using differen t abundance matching 
assumptions (e.g. Moste r et al.ll2013l) or different measured 
SMFs, leads to consistent results for the mass growth within 
~ 0.25 dex. The figure confirms that massive galaxies grow 
by a factor of 2 in stellar mass from z ~ 2 and a factor of ~ 5 
from z ~ 3, so that the typical stellar mass of the progenitors 
of ~ 10 1L2 M*/Mq galaxies is 10 10,5 at z ~ 3 and 10 100 at 
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z ~ 4. This mass growth track includes mergers, which occur 
at a rate of ~ 1.2 major (1:4) mergers/galaxy between z ~ 3 
and z ~ 0.5 (see section [731 for a detailed discussion on the 
effect of mergers). Since our sample is mass complete down 
to 10.0 from z ~ 3, a selection along the progenitors should 
not be affected by incompleteness. As described in lPatel et all 
dMl, we select galaxies along the growth track by picking 
galaxies in a given redshift bin within a narrow stellar mass 
bin of 0.3 dex around the corresponding mass for that red- 
shift. As also discussed in lPapovich et al ( 2014 ). this stellar 
mass bin is a reasonable trade off to account for the differ¬ 
ent predictions of different methods/SMFs and the scatter in 
the mass accretion histories. Table [T] summarizes the main 
properties of the selected sample at different redshifts. The 
redshift bins are selected to keep a comparable number of ob¬ 
jects in each bin (~ 400, except for the first and last bins) and 
as a tradeoff between comoving volume (~ 3.10 5 Mpc') and 
lookback time (0.5 — 1 Gyr). 

4. MORPHOLOGICAL EVOLUTION 

Figure [3] shows the evolution of the relative abundance of 
the different morphological types defined in § 12.21 selected 
along the mass growth track from figure [2] in 0.3 dex bins. 
The plot confirms the strong morphological evolution expe¬ 
rienced by the population of massive galaxies between z ~ 3 
and z ~ 1 essentially (see also e.g. Ivan Dokkum et al1l20111 
[Bruce et al .] l20H iMortlock et all 1201 3l> . About 60 - 80% 
of the progenitors of massive galaxies at z ~ 3 were irregu¬ 
lar disks (~ 40 — 50%) and mergers/irregulars (~ 20 — 30%), 
while the population at z < 1.0 is made at 80% — 90% by 
pure spheroids and galaxies with a classic bulge+disk struc¬ 
ture. Below z ~ 1, the well-known massive end of the Hub¬ 
ble sequence seems to be in place in terms of morphological 
mixing. Figure [4] illustrates this morphological transforma¬ 
tion with some example color stamps. Hence, considering all 
the progenitors of massive galaxies as an homogeneous fam¬ 
ily of objects when trying to infer their structural evolution, 
necessarily ignores the striking diversity of morphologies and 
the effect of morphological transformations. 

The inspection of the evolution for each morphology indi¬ 
vidually reveals some interesting trends. The fraction of pure 
spheroids is in fact roughly constant with redshift and repre¬ 
sents about ~ 30% of the population of massive galaxies at all 
epochs (only a slight decreasing trend is observed at z > 2.5). 
Most of the evolution is observed in the bulge+disk and the ir¬ 
regular disks populations which present more or less symmet¬ 
ric trends as clearly shown in the bottom panel of figured The 
latter goes from ~ 60% of the population of massive galaxies 
at z ~ 3 to roughly ~ 5% at z ~ 0.2. This decrease is mir¬ 
rored by the increase of the disk and bulge+disk populations 
which are almost inexistent at z = 3 and represent 50% of the 
galaxy population at low redshift. These trends suggest that 
most of the morphological transformations going on in the 
progenitors of massive galaxies go into one single direction, 
i.e from irregular/clumpy disks to more regular bulge+disk 
galaxies while the population of pure spheroids remains un¬ 
altered from z~2.5 and might follow an independent evolu¬ 
tionary track. 

The result might be an indication of two independent chan¬ 
nels for bulge growth in massive galaxies acting at very differ¬ 
ent timescales. Around ~ 30% of the population of massive 
galaxies at z ~ 0 was already made of bulges at z ~ 2.5 with 
probably an early ( monolithic ) fast collapse. The other half, 
however, have clearly a disk component and seem to appear 


gradually from z ~ 3 and z ~ 1 (~ 3 Gyr) through the mor¬ 
phological transformation of clumpy-irregular disks possibly 
through the migration of clumps and stabilization of the disks 
(e.g. Bournaud et al.ll20T4h . 

5. STAR FORMATION 

We now explore how the stellar populations, star formation 
rates and gas fractions evolve for each morphological type. 
Figures[5]to[6]show the evolution of the UVJ planes for differ¬ 
ent morphological types. Objects with different morphologies 
clearly populate different regions of the color-color plane as 
expected. Disk dominated galaxies (disks and irregular disks) 
tend to populate the star-forming region at all redshifts while 
pure spheroids are more concentrated towards the quiescent 
zone. Disk+spheroid galaxies lie between both regions. This 
confirms, that while a separation between passiv e and star- 
forming galaxies, as for example done by Ivan der Wei et all 
d20l4h . is clearly correlated with the morphology, it will not 
result in a clean separation of the morphological types and 
will mix bulges and disks. This is better seen in the left panel 
of figure |7] which shows the quiescent fraction for different 
morphologies, where quiescent galaxies are selected using the 
UVJ plane (red box in figures [5]to[6|. The average popula¬ 
tion is clearly quenched between z ~ 3 and z ~ 0.5 with the 
quiescent fraction rising from ~ 20% at z ~ 3 to ~ 8 0% a t 
Z ~ 0, in agreement with the findings of iPatel et alJ (1201 3l) 
and Panovi ch et all (120141. However, 90% of the disks and ir¬ 
regular disks are star-forming at all redshifts not showing any 
significant increase in the number of passive galaxies. On a 
similar vein, bulge+disk galaxies have a rather constant qui¬ 
escent fraction at all redshifts, close to 60%. The spheroid 
population however shows a clear increase going from a pas¬ 
sive fraction of 60% at z ~ 3 to almost 90% at z ~ 0. Given 
that the number density of spheroids remains roughly constant 
in the redshift range probed, this trend can be easily inter¬ 
preted as the same galaxies being quenched (within the limits 
of the abundance matching based selection). The increasing 
quiescent fraction observed for the overall population could 
then be explained as a combination of morphological transfor¬ 
mations of disk-irregular galaxies becoming disk+spheroids 
(as suggested by figure 0 and spheroids being individually 
quenched. 

A similar conclusion arises from figure [8] in which we 
plot the median star-formation rate (SFR) and specific 
star-formation rate (sSFR) for all morphologies. Different 
morphologies form stars at very different rates at all epochs, 
ranging from several hundreds of solar masses per year for 
the irregular and irregular disks to a few tens for spheroids. 
Generally speaking, objects with a significant bulge compo¬ 
nent tend to lie below the star formation main sequence at 
all redshifts (Whitaker et al. 2012 shown with stars in fig.0. 
The overall trend (black line in figure 0 is however a clear 
decrease of the SFR irrespective of the morphological type as 
predicted by sev eral models (e.g. iDekel & Birnboim; 120081 : 
iPeng et afll2010h and in agree ment with the evoluti on of the 
star formation main sequence (iWhitaker et alJl20l2h . 

Spheroids had a modest (compared to the average main se¬ 
quence at that epoch, i.e. I Whitaker et all 12012 1 star forma¬ 
tion activity already at z ~ 2 — 3 ( SFR ~ 50 M Q .yr~ l ) sug¬ 
gesting again that the formation of their stellar content oc¬ 
curred at earlier epochs and that they are in the process of 
quenching, i.e their star formation rate at z ~ 0.5 is almost 
0. We do observe however a significant increase of the aver- 
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Redshift 

Fig. 2.— Predicted mass growth of the progenitors of log(M * /Mq) = 11.2 galaxies from z ~ 4 from the lBehroozi et al] (20131 ) model. Error has show the 
errors on the median mass at a given redshift. 
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TABLE 1 

Summary of selected objects. (1) redshift range, (2) comoving volume probed in the corresponding redshift range considering 
THE CANDELS AREA, (3) LOOKBACK TIME INTERVAL, (4) STELLAR MASS RANGE, (5) TOTAL NUMBER OF OBJECTS, (6) NUMBER OF SPHEROIDS, (7) 
NUMBER OF DISKS, (8) NUMBER OF DISK+SPHEROIDS, (9) NUMBER OF IRREGULAR DISKS, (10) NUMBER OF IRREGULARS, (11) REMAINING GALAXIES 

WHICH INCLUDE UNCLASSIFIED, DISK+IRR+SPHEROIDS AND SPH+IRR 


age sSFR and above z ~ 2 it becomes larger than the thresh¬ 
old used by iBarro et all (120131) to define quiescent galaxies 
(log(sSFR\Gyr~ 1 ]) = —0.5). This increase is also accompa¬ 
nied by an increase of the scatter as also shown in figure[8] At 
z > 2, a significant fraction of spheroids are therefore actively 
forming stars, at similar rates than z ~ 1 main sequence disks 
(see also IBarro et al.ll2013h . 

Clumpy disks have rather high SFRs (> lOOMg.yr -1 ) at 
all epochs in which they are still abundant (z > 1 — 1.5) as 
well as disk dominated galaxies, roughly lying in the main 
sequence. Bulge+disk galaxies from roughly 50 Mg.yr” 1 de¬ 
parting from the star-formation main sequence. This suggests 
that, while the transition from irregular disks to disk domi¬ 
nated systems appears to be smooth without a big impact on 
the star formation activity, the morphological transformation 
between irregular and bulge+disk (i.e. the emergence of the 
bulge component) has to be accompanied by a decrease of 
their star formation activity and a departure from the main 
sequence. In other words, the emergence of the bulge and 


the stabilization of the disk in these objects tends to decrease 
significantly their SFR. This behavior is in line with the pre¬ 
dictions of several numerical simulations (e.g. iMartig et al.) 
1200 ^ ) which predict that the growth of a bulge in a turbulent 
disk can be sufficient to stabilize the gas disk and quench star 
formation (morphological quenching ). Another possibility is 
that the quenching which seems to follow the growth of the 
bulge is driven by the effe cts of a super massive black hole 
in the growing bulge (e.g. ISilk & Reesl Il998t iGranato et alJ 

l200l . 

To follow up on this idea, we look at the gas fractions 
through the existing correlation between the surface density of 
the star formati on rate and the cold gas through the Kennicut- 
Schmidt law (|Schmidtlll959t iKennicuttll 1 998l). As done in 
IConselice et ak (1201 3h and lPapovich et al.1 (120141) we use the 
following relation to infer gas masses: 

M gas = ( SFR \ 5/1 / R e \ 4/7 
6.8 x 10 8 Mg \lMQ.yr~ 1 J \lkpcJ 


(2) 
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FIG. 3. — Evo lution of the relative abundance of different morphological types as labeled between z ~ 0 and z ~ 4 along the mas s g rowth tracks 
from IBehroozi et all 120131) (see text for details). The shaded regions indicate the lcr error on the fractions co mpu ted following IGehrelsI 1198a) (see Section 
3 for binomial statistics; see also lMei et alJ|20()9|) . The top panel shows all the morphological types defined in !2.2l In the bottom panel, all the irregulars are 
combined in one class and all the disky galaxies in another one. 


which is then used to estimate the gas fraction as: 


f _ Mgas 

has ~ M gas +M t 


(3) 


Figure [9] shows the evolution of the inferred gas fractions 
for different morphologies. The average gas fraction de¬ 
creases monotonically with redshift from a value of ~ 60% 
t o ~ 20% at z ~ 0.5 as already reported in IPapovich et al.l 
(12014 ') for a slightly less massive sample. The trends differ 
significantly for different morphologies though. Spheroids, 
tend to have low gas fraction (~ 10%) at least from z ~ 2 


while, irregulars and disk-irregulars keep high gas mass frac¬ 
tions (50 — 60%). 

The increase of the gas fraction of spheroids at z > 2.5, even 
with the large uncertainties, is line with the idea of these ob¬ 
jects rapidly assembling at these epochs and consuming their 
gas reservoirs. On the other hand, the decrease of the gas 
content in disk galaxies, is again tightly linked with the emer¬ 
gence of the bulge component. While the decrease is rather 
smooth when no significant bulge is built, it becomes more 
dramatic for galaxies with a more predominant bulge (de¬ 
creasing from ~ 40% to ~ 20%). 
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Redshift 

Fig. 4.— Color stamps illustrating the evolution of the relative abundance of each morphological type along the main progenitors. Each column shows a 
different redshift bin, 0<z<l,l<z<2 and 2 < z < 3 from lest to right. The stamps are roughly built in the same rest-frame color using /814, /105, /125 
or /160 depending on the considered redshift. All stamps are normalized to the maximum pixel value. 
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Fig. 5.— Rest-fra me UVJ plane f or sph eroids (top panels) and disk+sph (bottom panels) at different redshifts as labelled. The red lines indicate the quiescent 
region as defined by I Whitaker et alj ( 120121) and gray points are all galaxies in the corresponding redshift/mass bin. 
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Fig. 6.— Rest-fra me UVJ plane f or disk s (top panels) and disk irregulars (bottom panels) at different redshifts as labelled. The red lines indicate the quiescent 
region as defined by I Whitaker et alj ( 120121) and gray points are all galaxies in the corresponding redshift/mass bin. 
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Fig. 7.— Quiescent fraction as a function of redshfit for different morphological types as labelled. In the left panel, the quiescent fraction is computed using 
the UVJ plane while in the right panel a threshold in sSFR (log(sSFR) < —1.5 [Gyr ~ 1 ]) is used. The trends are the same, but the absolute number of objects 
considered as passive changes depending on the definition used. 
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Fig. 8.— Median SFR (top) and sSFR (bottom) for different morphological types. Red: spheroids, yellow: disk+sph, blue:disks, violet:disk+irr, green:irr. The 
shaded regions show the 1 a uncertainty estimated through boo tstrapping. Black stars in the top panel show the position of main sequen ce galaxies at a given 
redshift according to the measurement of I Whitaker et alN2012f) . The red dashed-dotted line in the bottom panel shows the limit used bv lBarro et alj 420131) to 
define star-forming and quiescent galaxies. We only plot morphologies which represent at least 10% of the total population at a given epoch. 


6 . STRUCTURE 

We now move on to the study of the evolution of the 
structural properties of the different morphological types. 
Figure m shows the evolution of the effective radii, 
Sersic indices and axis-ratios. There is an average size 
increase of a factor of ~ 2 from z ~ 3, as a l ready pointed 
out b y many works (e .g iDaddi et all 120051: iTruiillo et all 
2006t iBui trago et al. 2008 1: Ivan Dokkum et all 120081: 
Damianovetal. 201 if |Cimatti et all 1201 2t iNewmanetalJ 
201 2i iHuertas-Co mpanv et alJl2013h . We do clear ly find two 
regimes in the size growth as also discussed by IPatel et all 
d2013l) for a similar selection. From z ~ 3 to z ~ 1.5, 
the average size of the whole population remains roughly 
constant and starts a sharp increase from z ~ 1.5 to z ~ 0. 
Recall that this differers from other works selected at fixed 
stellar mass (e.g. iNewman et alJl2012h because the selection 
is different. Adding the information of the morphological 
evolution discussed in §12.21 these two phases in the structural 
evolution are better explained. From z~3toz~1.5 there is 


a rapid morphological transformation of irregular disks into 
bulge+disk systems. Even though irregular disks are rapidly 
increasing their effective radii, their number density is also 
decreasing fast to be transformed into bulge+disk galaxies, 
which results in a decrease of the effective-radius because of 
the mass going into the central parts of the galaxy to build 
the bulge. As a result, both effects compensate to produce a 
flat size evolution. From z ~ 1 . 5 , the morphological mixing 
remains roughly constant and the average growth reflects 
simply the growth of the different morphological types. 
Interestingly, all dominant morphologies (spheroids, disks 
and disk+spheroids) at these redshifts do grow but the growth 
rate is different. While spheroids increase their effective radii 
by a factor of ~ 3 (~ 5 from z ~ 3), disks and disk-spheroids 
grow only by a factor of ~ 1 . 5 . The latter is roughly con¬ 
sistent with the expected growth of disks in galaxy haloes, 
i.e R e oc Z/(z) —1 (black dashed-dotted lines in fig. flOt which 
comes from the theoretical assumption that disks are formed 
with a fixed fraction of mass and angular momentum o f 
the parent halo (e.g. iMo et all 1 19981 iFerguson et al . 1120041 ) . 
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Redshift 

FlO. 9.— Median gas fractions or different morphologies. The black solid line show the average population while the different colors show different 
morphologies: red:spheroids, yellow, disk+sph, blue:disks. violet: disk irregulars and green irregulars. The shaded regions show the 1 — cr errors on the mean 
estimated through bootstrapping. We only plot morphologies which represent at least 10% of the total population at a given epoch. 


Spheroids grow at a faster rate as already pointed out by 
iHuertas-Companv et al.l (1201 3l) with a different selection, 
suggesting that some other mechanism takes place in these 
systems. 

Sersic indices also do increase on average from n ~ 1.5 
to n ~ 4, but again with different behaviors depending on 
morphology. The spheroids have n > 3 and they increase 
up to n ~ 5, confirming their bulge dominated morpholo¬ 
gies at all epochs. On the other side, irregular disks have 
very low n values (n ~ 1) while disk and disk+spheroids have 
rather constant intermediate values with (n ~ 1.5 for disks and 
n ~ 2.5 — 3 for disk+sph). This also confirms the validity of 
our morphological classification. Considering all these trends, 
the average observed increase of the Sersic index (black line 
in figure 1 1 (It is again better explained as a combination of 
morphological transformations from clumpy disks to regular 
systems which produces a growth of the bulge and an increase 
of the Sersic index together with the individual increase of 
spheroids. The general increase of the Sersic index is also ob¬ 
served by lShankar et al.l (i 2015) and might help to explain part 
of the evolution in the lensing profile of early-type galaxies. 

Axis ratios show little evolution with redshift but the abso¬ 
lute values change significantly with the morphological type. 
Spheroids have b/a values close to 0.8 while all the other mor¬ 
phological classes present values of 0.5 — 0.6, which again 
suggests that there are two families of objects following dif¬ 
ferent evolutions. The measured values are also in good 
agreement with measurements in the local universe for sim¬ 
ilar morphologies (iBernardi et al.ll20l3l) . 

In figure [TT] we analyze the total mass density profiles for 
different morphologies. As also done in iPatel et all ( 2013 ). 
we compute the median mass density profile using the best fit 
Sersic models at different redshifts and convert them to stel¬ 
lar masses by normalizing by the stellar mass of each galaxy 
(see lShankar et al.ll201 3l for details). This procedure is clearly 
a first order approximation since it neglects any gradient in the 
stellar populations which could definitely modify the shape of 


the profiles (specially for the star forming galaxies). The fig¬ 
ure clearly shows that spheroids are rapidly increasing their 
size at a faster rate than the average (bottom panel) with most 
of the action happening towards the galaxy outskirts through 
the addition of material. The central parts remain unaltered 
from z ~ 3 (changes in the inner 1 Kpc would not be detected 
given the PSF size). The evolution of the mass density pro¬ 
file for disks and disk+spheroids is less dramatic, resulting in 
a milder increase of their size, but the changes happen also 
at radii larger than 3—4 kpc. Interestingly, the clumpy-disks 
do show a significant evolution of their profile which could 
be interpreted as gas accretion happening in these objects and 
bringing material to the outskirts. These trends should how¬ 
ever be taken with caution, specially for the star forming pop¬ 
ulation, since firstly we neglect any stellar population gradient 
by construction when building the stacked mass profiles and 
secondly the profiles are obtained through single component 
fits which might not be well adapted to reproduce the irregular 
light distribution of clumpy galaxies. 

7. DISCUSSION: TWO CHANNELS OF BULGE 
GROWTH 

The results presented in the previous sections seem to point 
out two different evolutionary tracks for massive galaxies 
(/og(M*/M 0 ) ~ 11.2 ±0.15) and for the growth of their 
bulge component which are summarized in figure [12] As 
a matter of fact, the detailed analysis of the morphological 
properties of the progenitors from z ~ 3 shows that there are 
2 distinct families of galaxies with also different physical 
properties. 

7.1. The nugget track - Fast assembly 

About ~ 30% of massive galaxies had a spheroid morphol¬ 
ogy at z ~ 2.5 - with no disk component- and this fraction 
does not evolve down to z ~ 0. The quiescent fraction for the 
spheroid population is also rather high (~ 50%) and their gas 
fraction low (~ 10 — 15%) from z ~ 2 and increases between 
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Fig. 10.— Evolution of the median effective radius (left panel), Sersic index (middle panel) and axis ratio (right panel) for different morphologies. The 
red, yellow, blue, magenta and green lines show spheroids, bulge+disks, disks, irregular disks and irregulars respectively. The black line shows all galaxies 
irrespectively of their morphology. Error bars are estimated through bootstrapping. The squares show the values at z ~ 0 from IBemardi et al] 120141) and the 
dashed-dotted black lines in the left panel show the relation H(z)~ 1 normalized at the spheroids and disks z ~ 0.8 value. 
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Fig. 11.— Median stellar mass density profiles for different morphologies in different redshift bins as labelled. The top left panel, top right panel, middle left 
panel and middle right panel show the profiles for spheroids, disks, disk+spheroids and irregular disks respectively. The bottom panel shows the average profiles 
for all galaxies. 
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z ~ 2 — 3 which suggests that these galaxies are rapidly assem¬ 
bling at z ~ 2.5 and above. The fraction of passive spheroids 
rises to almost 100% at z < 0.5 indicating that they are in 
the process of quenching in the epoch probed by this work 
but without significant alterati ons of their morphologie s, in 
agreemen t with the findings of iBarro et alJ (120131120141) and 
iMei et all (120141) for dense regions. Their SFRs remain in¬ 
deed well below the main-sequence of star formation at all 
epochs (fig. A3. The fast quenching is also accompanied by 
a rapid growth of their effective radii by a factor of ~ 5, com¬ 
pared to a factor of 2 growth in stellar mass, with most of the 
action happening in the galaxy outskirts (R > AKpc). The Ser- 
sic index also increases from n ~ 3 to n ~ 5, even though it 
remains rather high at all epochs, confirming their bulge dom¬ 
inated morphology. They are th erefore very simila r to th e 
dense-core galaxies identified by Ivan Dokkum et all (120141) . 
The increase of the Sersic index, is however not coupled to the 
gas content, at variance with what is observed for the average 
population (iPanovich et al.ll2014h which suggests an external 
driver. The fact that the number density and the morpholo¬ 
gies do not change supports the idea that the reported growth 
is indeed an individual growth of these objects and that it is 
not driven by the morphological transformations or quench¬ 
ing of new galaxies (progenitor bias). Otherwise, we would 
expect an increase of their abundance, since it is very unlikely 
that these galaxies will transform into another morphological 
class (although this cannot be fully excluded as discussed in 
section E2>. An independent test for this statement would 
be a detailed analysis of the stellar population ages for this 
particular population which will be addressed in forthcom¬ 
ing work with higher resolution SEDs. Recall, that the size 
growth is even larger than the one measured for the overall 
population (i.e. factor of 2-3) which is in fact a convolution 
of different mechanisms as described below. The properties 
of these galaxies are therefore consistent with a formation of 
the bulk of their stellar populations at high redshifts through 
violent disk instabilities or gas rich mergers (at z ~ 2.5 their 
gas fraction is high and also their sSFR) which create d their 
dense cores (IBarro et alJl2014t iDekel & Burkertir2014l) and a 
subsequent growth by the addition of material in the outskirts 
through, possibly, minor mergers. Figure [12] summarizes the 
inferred evolution of these objects in the mass-size and M, - 
SFR planes from z ~ 3. It is worth emphasizing, that the size 
growth for these objects at later epochs (z < 1) is still very 
pronounced, i.e. a factor of ~ 2 with a minimum stellar mass 
growth. The global trend is nevertheless still compatible, at 
first order, with a minor merger driven growth as predicted 
by numerical and semi-empirical models (e.g. Shankar et alJ 
120151 for slightly larger stellar masses). A more detailed com¬ 
parison with the models predictions is however required at 
this stage. 

7.2. The clumpy track - Slow assembly 

The remaining ~ 60% of the population is made of irreg¬ 
ular/clumpy disks at z ~ 3 which experience a rapid mor¬ 
phological transformation between z = 3 and z = 1 to give 
birth to very massive disks with a small bulge (~ 20%) and 
to 40% of galaxies with both a prominent bulge and a disk 
component. The evolution of the effective radii is more mod¬ 
erate than for spheroids and scales roughly with //(z) -1 , the 
expected growth of disks in DM haloes. Figure [12] summa¬ 
rizes the inferred evolution of these objects in the mass-size 
and M*-SFR planes from z ~ 3. The transition between the 
clumpy-irregular morphologies to more Hubble sequence like 


galaxies happens mostly at z > 1 ■ Clumpy-irregular disks are 
characterized by high SFRs (> lOOMg.yr - ', slightly above 
the main-sequence at that epoch), high gas-fractions (~ 60%) 
and low Sersic indices (n ~ 1). Some of these objects (~ 1 /3) 
will experience a smooth transition to become massive spirals 
with low B/T fractions. As a matter of fact, the properties of 
both families are very similar, in terms of gas fractions, Ser¬ 
sic index, SFRs and also effective radii. The other ~ 2/3 will 
build a more prominent bulge (n > 2.5) which roughly corre¬ 
spond to B/T of 50%-75% (iBruce et al.ll2012l) . The building- 
up of a larger bulge component results in a decrease of the ef¬ 
fective radius following the concentration of mass towards the 
central regions. The emergence of the bulge is also tightly cor¬ 
related with the decrease of the star formation activity and the 
decrease of the gas fractions which go from 50% to ~ 10%, 
which make them depart from the star-formation main se¬ 
quence tfig.n~2l>. This evolution is consistent with the predic- 
tions of several n umerical simulations (e.g. Bournaudl l2015l: 
Marti g et al.l l2009h which show how the growth of a bulge 
through clump migration is followed by a decrease of the 
star-formation activity (morphological quenching) although 
the effect of feedback from a SMBH in the growing bulge 
could also produce similar effects. As a matter of fact, AGN 
feedback is k nown to contribute to the quenching of s tar for¬ 
mation (e.g., ISilk & Reesl 119981 : iGranato et all 12004]) and it 
is also known to correlate with the mass of the bulge (e.g., 
iKormendv & Hol2013l ) and therefore could also help explain¬ 
ing the quenching of bulge+disk systems which seems to be 
associated with the growth of the bulge. 

7.3. Major mergers 

This 2 track-scenario is obviously not the only possible ex¬ 
planation for the trends we observe. Namely, the constant 
number of spheroids could also be a result of clumpy gas rich 
galaxies being transformed into spheroids following a major 
mer ger event and spheroids regrowing a disk at the same rate 
(e.g. iHammer et al]|2009t IHopkins et al.ll2009l) . Since the in¬ 
ferred gas fractions of irregular systems is high (~ 40%), this 
is a plausible option. The scenario requires however a fine- 
tuning to keep these two effects (formation of spheroids and 
disk regrowth) at the same rate and also a high (major) merger 
fraction to keep producing spheroids. There have been sev¬ 
eral measurements of the major merger rate (1:4) of massive 
galaxies (M*/M 0 ~ 10 11 ) in the recent literature. From the 
observational point of view, iLotz et akl (1201 ll) measure 1.6 
mergers/galaxy between z ~ 0 — 3 (extrapolating the quoted 
redshift evolution). This is in rather good agreement with 
iMan et akl(l2012l) who find Lj_ mergers/galaxy i n the same pe¬ 
riod and also with iLopez-Saniuan et aD (120121) . [Bluck et al. 
(120121120091) find a larger fraction (1.7 mergers/galaxy only 
between z = 1.7 and z = 3 ). Abundance matching based 
measurements (e.g. IHopkins et al.l I 20 I 0 I) also find similar 
numbers (1.7 mergers/galaxy between z = 0 and z = 3), just 
as SAMs (e.g . ICon selice et akl 1201 4l) and numerical simu¬ 
lations (e.g. iKavirai et all 1201 4lk In this work we use the 
iBehroozi et al.l d2013l) model, which predicts 1.2 major merg¬ 
ers/galaxy along the mass growth track shown in figure |2] 

Considering these different measurements, it is certainly 
safe to assume that, on average, each galaxy in our sample ex¬ 
periences a major merger event in the redshift range explored. 
Assuming then that each merger event is enough to change 
the morphology, it is indeed possible to explain the number 
density decrease of irregular disks by mergers followed by 
disk rebuilding. We notice however that the IHopkins et al.l 
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d2009l) simulations focusing on disk rebuilding, predict that 
only equal mass mergers are able to create bulge dominated 
systems. Lower mass ration mergers (1:2-1:4) tend to cre¬ 
ate disk dominated systems which continue forming stars. In 
that respect, if mergers+disk rebuilding is the dominant chan¬ 
nel, we would expect at all epochs an increasing fraction of 
star-forming regular disks and a minor fraction of quiescent 
bulge-dominated systems. The opposite is actually observed. 

In addition to this, we do observe that the median ages of 
the spheroids in our sample estimated through SED fitting, 
monolithic ally increase from ^0.5 Gyrs at z > 2 to ~2.5 Gyr 
at z ~ 0.5 (roughly consistent with the time between these 2 
redshift bins). On the other hand, the ages of disks dominated 
systems tend to stay rather young (< 1 Gyr) at all epochs due 
to the sustained star formation. If the dominant process to cre¬ 
ate the 30% spheroids we observe is merging, we would not 
expect a strong increase of the ages of spheroids which we 
actually seem to observe. Given the known degeneracies af¬ 
fecting age determination from broad band photometry, these 
trends need to be taken with caution. However they point to¬ 
wards an early formation of the spheroid population. A fur¬ 
ther check of the proposed bulge growth tracks would there¬ 
fore imply accurate age estimations of the bulge components 
in the different galaxy types. This requires bulge-to-disk de¬ 
compositions and high resolution SED fitting of the different 
components which is on-going. 

8 . SUMMARY 

We have analyzed the morphological, structural and star- 
formation properties of the progenitors of massive galaxies 
(M. t /M ( .j ~ io 1L2±0 - 3 ) from z ~ 3. The progenitors are 
selected using abundance matching to take into account 
the expected mass growth in the redshift range probed in 
this work. The latter selection is a key point of the present 
work and is clearly subjected to important assumptions 
(i.e. halo mass functions, galaxy stellar mass functions and 
halo occupation distribution functions) as explained in the 
previous sections. It is worth emphasizing however that 
selecting galaxies at fixed stellar mass (i.e. assuming the 
extreme scenario in which galaxies do not grow in mass) 
results in very similar trends that the ones reported here. 

The main new ingredient of this work is the addition 
of accurate visual-like morphologies which helps better 
understanding the different evolutionary tracks leading to 
the present day Hubble sequence. We have defined 5 main 
morphological types which quantify the presence or not 
of a bulge/disk component and the presence or absence of 
irregularities in the light profile. We then have explored, the 
abundances, star formation rates, quiescent fractions, gas 
fractions and structural properties for each morphological 
type. 

Our main results are the following: 

• The morphologies of massive galaxies significantly 
change from z ~ 3. At z < 1, these galaxies are made by 
40% of pure spheroids, 40% bulge+disk galaxies (early 
spirals and lenticulars) and 20% of massive disks. At 
z ~ 3 there is still a 40% of spheroids, but the remain¬ 
ing 60% is made of irregular/clumpy disks or disturbed 
galaxies. Most of the morphological transformations 
take place at z > 1. 

• As reported in previous works, the overall population of 


massive galaxies rapidly quenches from z ~ 3 to z ~ 0, 
i.e. the quiescent fraction increases from ~ 20% to 
~ 80%, the median SFR decreases from ~ lOOMg.yr - 1 
to ~ 25Mg.yr -1 and the gas fractions go from ~ 40% 
to ~ 15%. When inspected at fixed visual morphol¬ 
ogy the trends are very different. The quiescent frac¬ 
tion in the spheroid population is already high at z ~ 3, 
i.e. 60% and increases to an almost 100% value. The 
quiescent fraction for disks and disk irregulars remains 
low (< 20%) at all epochs while for bulge+disk ob¬ 
jects the fraction appears to be constant too at a value 
of 40 — 50%. These trends suggest that the overall in¬ 
crease of the quenching fraction for the whole popula¬ 
tion can be explained by a combination of the quench¬ 
ing of the spheroid population with the morphologi¬ 
cal transformation from clumpy/irregular disk to early- 
spiral/SO. 

• When considering the overall population, without mor¬ 
phology distinction, we measure an increase of the av¬ 
erage effective radius by a factor of 2 — 3 as well as 
an increase of the Sersic index from n ~ 1.5 to n ~ 4, 
as reported in the recent literature for similar selec¬ 
tions. The evolution of the average mass density pro¬ 
file is also in agreement with an inside-out growth. The 
evolution of the average size seems to have 2 differ¬ 
ent regimes, from z ~ 3 to z ~ 1.5, there is almost 
no significant change of the effective radius while the 
bulk of the growth happens from z ~ 1.5 to z ~ 0. At 
fixed morphologies, spheroids do grow by a factor of 
5—6 from z ~ 3 and increase their Sersic index from 
n ~ 3 to n ~ 5. On the other hand, irregular disks and 
disks grow by a factor of ~ 1.5 and keep a rather con¬ 
stant Sersic index (n < 2 for disks and disk+irr) and 
n ~ 2.5 — 3 for disk+bulge galaxies. The two differ¬ 
ent phases in the average growth are better explained if 
morphological transformations are taken into account. 
In the first phase, there is a rapid morphological trans¬ 
formation from clumpy disks to bulge+disk galaxies 
which results in a slight decrease of the effective radius 
as a consequence of the mass concentration towards the 
inner regions of the galaxy. Even though, spheroids 
and clumpy disks increase their size in this period, the 
rapid decrease of the number density of the latter seems 
to compensate this growth and results in no evolution 
of the average size. During the second phase, from 
z ~ 1.5, the morphological mixing remains roughly 
constant, but the size growth increases by a factor of 
~ 2 — 3 on average. This growth is therefore better 
explained by the individual growth of disk/disk+bulge 
galaxies which grow by a factor of ~ 1.5 and the growth 
of the spheroids which increase their effective radius by 
a factor of ~ 4. 

The above results suggest to different channels for the bulge 
growth in the massive end of the Hubble sequence: 

1. A nugget track (fast assembly) followed by 30 — 40% 
of the population of massive galaxies. Galaxies formed 
that way, formed the bulk of the stars at z > 2.5 
and also acquired their spheroidal morphology at these 
early epochs possible through violent disk instabilities 
(and/or mergers) which rapidly bring gas into the cen¬ 
tral parts. At z < 2, they have already low gas fractions, 
low SFRs, high Sersic indices and ~ 60% of them are 




18 


Huertas-Company et al. 





3.0 


2.5 


A 2.0 

o: 

Ll. 

c/t^ 

cr> 

€ 1.5 


1.0 


0.5 



z ~ 2 

Z ~*z's 5 ;;.'.'.' 

z ~.2-0"” z ~i '5-.^ 

y ■' Bat. 1 z~1 


Z- 


\A-I 






10.6 10.8 11.0 11.2 11.4 

Log(M./M 0 ) 


Fig. 12.— Expected evolution in the mass-size (top) and M* -SFR (bottom) planes of the two channels of bulge growth (see text for details). The left panels 
show the evolution of spheroids. The right panels show t he ev olutionary track of clumpy disks. The red, yellow and blue solid lines in the top panels show 
the z ~ 0 median mass-size relation from lBernardi et all (2014) for ellipticals, early-spirals and late-spirals respectively. The dashed lines in the bottom panels 
indicate the median star-forming main sequence at different redshifts from lWhitaker et SI (2012 ) 


classified as quiescent. They are however very compact 
with median effective radii of ~ 5 kpc. Between z ~ 3 
and z ~ 0.5 they practically completely stop forming 
stars while they increase their size by a factor of ~ 5 
and their Sersic index from n ~ 3 to n ~ 5 keeping their 
global visual aspect unaffected. The growth is decou¬ 
pled of the gas content and the SFRs which remain low 
at all epochs and mostly happens in the galaxy outskirts, 
suggesting an ex-situ driven growth. 


2. A clumpy track (slow assembly), followed by 60 — 
70% of the population of massive galaxies at z ~ 
0. These galaxies were clumpy/irregular star-forming 
disks (SFR > lOOM 0 .yr _1 ) at z ~ 2 — 3. From z ~ 3 
to z ~ 1 they experience a rapid morphological trans¬ 
formation leading to relaxed systems (at least in terms 
of their visual aspect) and to the emergence of a bu¬ 
gle component of variable size (~ 2/3 seem to de¬ 
velop a large bulge component while the remaining 1 /3 
keep a disk dominated morphology). The morphologi- 
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cal transformation is accompanied by a decrease of the 
SFR (reaching ~ 50Mg.yr _ 1 ) and the gas fraction (go¬ 
ing down to ~ 15%), more dramatic for galaxies de¬ 
veloping a larger bugle as well as by an increase of 
the Sersic index (from n ~ 1 to n ~ 2.5 — 3) and a de¬ 
crease of the effective radius because of the mass be¬ 
ing concentrated towards the inner regions. This is in 
good agreement with the predictions of numerical sim¬ 
ulations in which the bulge component is built from the 
migration of clumps and the stabilization of the disk 
results in a decrease of the star formation rate ( mor¬ 
phological quenching) although the possible effect of a 
SMBH should also be considered. Major merger events 
followed by a disk rebuilding event could also con¬ 
tribute to transform irregular systems but it is unlikely 
to be the dominant channel in the mass range explored 
in this paper. Below z ~ 1, the well known massive 
end of the Hubble sequence is in place and the galaxy 
properties change only marginally. Their effective radii 
grow in fact at a rate roughly consistent with H(z)~ 1 
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